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3 Temperature dependence of the lower critical field H c i (T) for the superconducting magnesium 
diboride, MgB 2 , is studied in the vicinity of T c by using a two-band Ginzburg-Landau (G-L) theory. 

O 

Theoretically calculated Hci (T) near T c exhibits a negative curvature. The results are shown to be in 
a good agreement with the experimental data. In addition, two-band G-L theory calculations give a 

/t(T) 

temperature dependence of the Ginzburg-Landau parameter k(T) = , which varies little with 
temperature as in a similar manner of the microscopic single-band BSC theory. 
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I. INTRODUCTION 

The recent discovery of superconductivity in an intermetallic compound of magnesium 
diboride, MgB 2 , [1] is interesting for some reasons. First, MgB 2 holds the record of the highest 
critical temperature at around 40 K among the simple binary compounds of superconductors. The 
second, the material shows a pronounced isotope effect of boron, which suggests that MgB 2 is a 
phonon mediated superconductor [2]. In addition, the compound has been found particularly 
attractive because it is generally thought that MgB 2 has a potential to replace conventional 
superconducting materials in the electronic applications in terms of large critical current densities [3] 
and transparency between grain boundaries to the current flow [4]. MgB 2 is of a type II 
superconductor with the estimated values of a coherence length %(0K) , penetration depth X(0K) and 

G-L parameter, k(0) being 5.2 nm, 125-140 nm, k(0)=26, respectively [5,6]. 

Microscopic mechanisms of superconductivity in MgB 2 and the discussions of its 
superconducting properties implicitly involve the employment of Isotropic Single Band (ISB) model 
[7]. However, band structure calculations suggest that the origin of the superconductivity is based on 
the hole-like boron o band [8] and % band [9]. Hence, two-band Eliashberg theory [10] is generally 
assumed as complementary to the ISB model. Recently two-band Eliashberg model was applied for 
determining the behaviour of H c2 in MgB 2 [11]. 

The temperature dependence of the lower critical magnetic field H c i exhibits a negative 
curvature near critical temperature T c as reported experimentally in Refs. [12-14]. For a review, see 
also the recent work of Ref. [15]. Although there has been some extensive theoretical works on 
MgB 2 , to our best of knowledge, there is no systematic model to explain the temperature dependence 
of H c i, observed in MgB 2 . In a previous communication, we presented the temperature dependence 
of H c2 (T) by using the two-band G-L theory [16]. In this paper, we show that the presence of two 
order parameters for two bands in the theory yields a negative curvature in the vicinity of the critical 
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temperature for MgB 2 superconductor. The calculated temperature dependence of the lower critical 
field H cl (T) is in a good agreement with some experimental data. 



II- THEORY 



We write the G-L free energy functional for two-band superconductors with two coupled 
superconducting order parameters of the following form similarly as in [16,17], 

F SC [ < P 1 ,^ 2 ] = \d i r(F, +F n +F 2 +H 2 /&r) (1) 
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Here m i denotes the the effective mass of the carriers belonging to the band i (i=l,2). Fj is the free 
energy of the separate bands. Fi 2 is the interaction energy term between the bands. The coefficient a 
depends linearly on the temperature T as OCj = Yj (T -T c j ), while the coefficient (3j is independent of 
temperature. Here, Yj is the proportionality constant. The quantities £ and Ei describe the interband 
mixing of two order parameters and their gradients, respectively. H is the external magnetic field , 
3> () is the magnetic flux quantum and A is the vector potential. In Eqns. (2) and (3), it is assumed 

i i2 

that the order parameter *Pj | are slowly varying functions in space. 

For temperatures near T c and magnetic fields slightly larger than H c i, the influence of the 
field on the modulus of the order parameters *P, and m 2 can be neglected and we assume I^Uconst 



hP 2 | =const . Then the representing wave function v I , i .(r)can be written as 



*P,-(r) = | y ¥ i |expO'^(r)) . Here, (j) j (r) axe the phases of the order parameters and the G-L free 
energy functional (1), can be rewritten as 
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r h 2 dd)i 2tc - ? # 2 d<b 271 - ? 

— mCTX-p--— A) 2 +— n 2 (T)(4 2 --— A) 2 + 
8m! dr <£ 8m 2 dr <J> 

+ e( ni (T)n 2 (T)) 1/2 cos((h-(h) + 

+ 8 1 (n 1 (T)n 2 (T)) 1/2 cos# 1 -(^ 2 )(^--^A)(^ 2 -^A) + H 2 /87l 

dr $ dr *0 

V 



(4) 



Where, n x (T) = 2 and n 2 (r) = 2 *P 2 are the densities of superconducting electrons for the 
corresponding bands, respectively. The temperature dependencies of m (T), n 2 (T) and (0j -0 2 ) are 
defined by the equilibrium value of order parameters \ x ¥ 1 1 and |*P 2 1 , which satisfy the G-L equations 
(see Eqn. A2, A3). 

The equations determining the equilibrium values of magnetic field and phases of the order- 
parameters would be obtained by the rninimising the free energy functional (4) with respect of the 

vector potential A and the phases (j)j , (j) 2 . The equation for the vector potential takes the form 
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By using the relevant Maxwell equation, VxH = — J (for the magnetostatic case), Eq. (5) yields to 

c 

the London equation ( by taking into account the equilibrium value of the differences of phases 
0! - (j) 2 , see A7 and A8) of the form 



d H 
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where X is the London penetration depth of the following form 
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It is well known that, the lower critical field H c i can be obtained as in Ref. [18] 



In k(T) . 
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We then introduce a dimensionless lower critical fi eld h. = — - with 

H cl (0) 

H , (0) = — ^— — — ( — - — + — - — ) . We rewrite the normalised lower critical field as: 
c 2 Pimi P 2 m 2 



h cl =B(T)\nK(T), 
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B{T) = -—^—ris 2 + *(t - t c1 ) 2 + 2^r7£ 2 (t - t c1 )) 



e 2 + (2-T cl -T c2 )e 

^£ 2 D(t-t c2 ) + (t-t c1 ) 3 



(10) 



7 

el,e2 



Temperature dependence of the normalised G-L parameter k(T) in two-band SC is given as 



, X = 



Yimi ,0=T-l,Tl= £lTcY2m2 
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Here, the upper critical field h c2 of two-band superconductors was calculated as in [16] and can be 
given as 



h c2 =a- 1 {-e-c +(A6 2 +B0 + c o ) U2 ) 
The parameters in last equation are the same as in [16]. 



(13) 



III- RESULTS AND DISCUSSIONS 

The experimental data for the lower critical field hd of MgB 2 can be fitted well with high 
accuracy by the simple expression of the form 

*d=*e'l(-0) lH \ (14) 
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In figure 1, we plot the experimental data of Li et al. [12] with open circle symbol and 
Sharoni et al. [14] with filled circle symbol for H X (T) in MgB 2 versus the reduced temperature T/T c 
with H cl (0) = 290 Oe and 450 Oe respectively. We also show the theoretical fits to the data by 
using Eqns.(9) and (14). The best fit to the experimental data was obtained by using Eqn.(14) (shown 
by solid lines) with fitting parameters h* cl = 0.9 ,a = 0.23 . As shown in [16], the upper critical field 

can be fitted by the similar formula with x+a ,a = 0.23 In Fig.l, the dotted line displays the two- 
band G-L fitting by using Eqn.(9) with D=1.35. The temperature dependence of the two-band G-L 
k(T) 

parameter — — is shown in Fig. 2. In the calculations, G-L parameters for the fitting h c2 are the 
same as in [16] 

A = 0.59, B = -0-11, c = 0.20,x = 3,|l H c2 (0) = 15.55 T,T cl =20 K,T c2 =10K. 

k(T) 

Calculations show that the function has weak temperature dependence as in a similar way of 

K(0) 

the single-band BCS theory [18] implying that the temperature dependence is determined mainly by 
the temperature dependence of the penetration depth. On the other hand, as followed from the 
equation (7), h c i is determined dominantly by the small mass of the carriers in contrast to the upper 
critical field, h c2 . The term in Eqn. (7) for the intergradient interaction between the bands may also 
contribute to the temperature dependence of the penetration depth. Inclusion of the term with £j 
enhance the negative curvature near critical temperature (dashed line) in similar manner for upper 
critical field. Thus, in two-band G-L theory, the lower critical magnetic field exhibits negative 
curvature in contrast to ISB GL theory. Note that the single band G-L theory would only give a 
linear temperature dependence for H c i [18] only at temperatures close to T c . 

As shown by Eqn. (A4) , the critical temperature of around T c =40 K can be obtained from 

T c i=20 K and T c2 =10 K approximately with the interaction parameter 8 =3/8. As followed from the 
eq. (11) lower critical field h c i is governed by the parameter 8j, 8 . Our fitting corresponds to T| = - 
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0.16. We also take into account that the average Fermi velocity over hole tubular band is v F1 =4.5 
10 7 cm/s, while for the remaining last Fermi surface sheets we can choose v F2 =11. 10 7 cm/s . 

Note that the theoretical results of two band superconductivity is in a good agreement with 
the experimental works of [12] and [14]. However, there are some experimental works that are not 
well fitted to our results at high temperatures [19] and at low temperatures [13]. There may be 
several reasons for the inconsistency. First, Joshi et al. [19] has presented experimental 
measurements of H c i well below T c and found a linear dependence on temperature. It is generally 
accepted that H c i becomes vanishigly small near T c and its measurement is almost impossible at 
temperatures in the vicinity of Tc. We also note that our theoretical results for hc2 [16] are in 
agreement with their experimental data of Ref. [19]. In addition, the measurement of H c2 (T) in Ref. 
[13] also shows positive curvature near Tc. Their data for H c i(T) are in a reasonable agreement with 
the data of Sharoni et al. [14] at temperatures between 20 and 40 K, while their data differ at 
temperatures below 20 K. In our opinion, the discrepancy may result from the sample specifications, 
since optimisation of the sample preparation is still needed to make high quality samples. Two-band 
G-L theory gives better temperature dependence than single-band G-L theory. This may support the 
proposals that MgB 2 is the first two-band superconductor [8,9,11,20] 

In summary, we have shown that the negative curvature of the lower critical field can be 
explained to a reasonable extent by seeking analytical solutions to the two-band G-L theory within an 
approximation. An exact fitting between experiement and the theory would necessitate exact 
solutions to the theory. However, the theory itself is phonemologic and it can explain the temperature 
dependence to some extent only at temperatures close to T c . Further work is therefore needed to 
elucidate the nature of superconductivity and the interesting physical properties observed in MgB 2 at 
both microscopic level and the macroscopic level. 
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APPENDIX 

By minimization of the free energy of Eqn.(l) 

SF/9¥t=0 , 8Fld¥* 2 =0 (Al) 
we obtain the usual basic equation for the description of the two-band superconductivity. For 

simplicity, we assume a vector potential A = (0,Hx,0) in one dimension, then linearisation of 
Eqn.(4) yields 

' (TT - 4)^1 + «i (T)% + ^2 + fii (4t - TT^ + JW' = (A2) 



4m, dx / s " dx I 

H ~ <l ' +a 2 (T)V 2 +£,(4^-^-)*, + P 2 Vt = . (A3) 



4m 2 dx Z s dx / 

Where lg = (/zc/2eH) is the so called magnetic length. The coefficients 8i and e are determined by 

the microscopic nature of the interaction between different band electrons. Their signs can be taken 
arbitrarily. If the interband interaction is ignored, the Eqns. (A2) and (A3) are decoupled into two 
ordinary G-L equations with two different critical temperatures. In general, regardsless of the sign of 
8, the superconducting phase transition occurs at a well defined temperature exceeding both T c i 
and T c2 , which can determined by the equation 

a 1 (T c )a 2 (T c ) = e 2 . (A4) 
In the absence of any magnetic field, the equilibrium value of the order parameters are given as 

^ = g^CO^CO-g 2 ) 
e 2 j8 1 a 2 (r) + j8 2 a 1 3 (T) 

x V 2 = a li T ) fa{T)a 2 (T)-e 2 ) 
e 2 p ia2 (T) + l5 2 af(T) 

Phases in the equilibrium state satisfy the following condi tions 

cos(0 t -0 2 ) = 1 or 1 -0 2 =2tcti , if £<0 (A7) 

cos(0! - 2 ) = -1 or 0, - 2 = (2n + l)n , if e >0. (A8) 
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Figure Captions 



Figure 1. Temperature dependence of the upper critical field H c i(T). Open circle symbols show the 
experimental data of Li et al. in Ref. [12], while open square symbols show the experimental data of 
Sharoni et al. in Ref. [14]. The solid line shows the calculations from Eq.(14). Filled circle symbols 
with dotted line represent the two band G-L theory (8j =0). Triangle symbols with dashed lines are 
the two-band theoretical calculations of G-L theory ( 8j = 0.0976 ) involving interaction between 
bands. 

Figure 2. Temperature dependence of Ginzburg-Landau parameters K . The solid line represents 
calculations from the theory with 8j = and the dashed line is for 8j = 0.0976 . 
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